Terahertz wave emission from mesoscopic crystals of BSCCO 
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We study the Terahertz (THz) emission mechanism in BSCCO with an external magnetic field 
theoretically. The total power is proportional to the square of the number of layers and the fre- 
quency. Under the small size limitation, the angular distribution is similar to that in the dipole 
emission. The present theoretical results for THz power output are in the same order of magnitude 
as the experimental ones, in contrast with previous simulations. A series of size-dependent intrin- 
sic frequencies are also observed. When the external frequency induced by the external current 
resonates with intrinsic ones, the distinct THz emission is observed. 

PACS numbers: 74.50.+r, 74.25.Gz, 85.25.Cp 



Highly anisotropic high-Tc superconductor 
Bi2Sr2CaCu2 Os (BSCCO) is considered as a stack 
of intrinsic Josephson j unctions (IJJs) on atomic scale 
p], [2J. This materials may provide a new way to gen- 
erate terahertz (THz) electromagnetic (EM) wave and 
therefore possibly fills the "THz gap" @, |, i, !, 0, B 0, 
one area of the EM spectrum between microwave and 
infrared frequencies, which has not yet been exploited. 
The THz EM-wave emission based on AC Josephson 
effect of IJJs has been detected by Ozyuzer et al. [J], 
and also studied extensively in the theoretical side. 
Some mechanisms of the THz emission have been 
proposed [5|, l3| recently. However, these mechanisms can 
not bridge the big gap for the power in simulations 
and experiments [4|, 

In this paper, we successfully develop a new boundary 
conditions, which explains the dipole angular distribu- 
tion and predicts that the power is proportional to the 
square of layer number and frequency. By solving the 
sine- Gordon equations numerically with a proper bound- 
ary conditions, we observe a stable power output in some 
parameter regime. The present power density is in the 
same order of magnitude as the experimental results. Fi- 
nally, the detail mechanism is suggested. 

The plasma frequency of single Josephson junction is in 
subterahertz wave band. There is a collective cavity res- 
onance mode in the stack of IJJs, which is crucial to en- 
hance the power and to increase the frequency from sub- 
THz to THz. IJJ is formed naturally in BSCCO where 
Bi-Sr-0 layers between the superconducting Cu02 layers 
acts as a non-conducting barrier. It can be described by 
the well known coupled sine-Gordon equation [1] 
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where x' = x/X c and y' = x/X c are the dimensionless 
positions in the a-b plane, t' 
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FIG. 1: (Color online). Schematic view of BSCCO source 
(left) and the equivalent radiation sources (right). 



insulating layer thickness. (3 = 47rcr c A c /cy^ is the resis- 
tive dissipation factor, where e c and a c are the dielectric 
constant and the conductivity of the insulating layers, re- 
spectively. The abbreviation operator of A^ 2 ^ is defined 
as A( 2 )/(Z) = /(Z + l)-2/(Z)+/(Z-l). 

It is known that the boundary conditions are very cru- 
cial for the wave emission. The BSCCO source and out- 
side media region are separated by a boundary surface. 
On the interface, the surface integrals of E and H are 
constant, except in the red region shown in the Fig. [H It 
is an effective radiation section where H parallel to y-axis 
and E to x-axis are determined through <p t — 2eDE z /h 
and (j) x = 27rjj J oDHy/^o respectively. So the EM-wave 
is described by Maxwell equations outside the BSCCO 
region and sine-Gordon equations inside the BSCCO re- 
gion. To evaluate the far field, the original problem could 
be transformed to an equivalent one described by Love's 
Formulation [Toj] . In this way, the original IJJ region is 
replaced by a free space, the BSCCO source is removed, 
and equivalent electric and magnetic currents given be- 
low are suspended on the boundary surface, 

J = n x H, 



J n 



-n x E, 



where n is an unit vector normal to the effective radiation 
section, J and J m are the equivalent sources producing 
the same E and H in the region outside. According to 
the superposition principle of the electromagnetic waves, 



2 






(b) 

V / 


"v. \ 


},. 

}. O 


O 










J 0, 4 









Sx 



FIG. 2: (Color online), (a) Angular distribution of the radi- 
ation energy in ab-plane. The present analytical results are 
denoted by the solid curve and the previous simulation ones 
by the open circles, (b) Angular distribution of the radiation 
energy in ac-plane. 



the E and H in far field are the sum of contributions from 
electric sources and magnetic sources respectively. The 
angular distribution of total power is determined by the 
radiation from both J and J m . For the electric sources, 
the vector potential in outside region is generated by J, 



jft \ e ik\r-r'\-ujt 
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The average energy flow density can be evaluated by the 
time averaged Poynting vector 



S{r)= 32.2c /<L ^ fln 



(3) 



It is Eq.(j3]) that can account for a lot of experimental 
and numerical observations, such as the angular distribu- 
tion of the emission, total power which is proportional to 
iVj, and considerable power discrepancies between simu- 
lations and experiments. 

The analytical results for the angular distribution of 
the radiation energy by Eq. (3) is plotted in Fig. [2] 
with the solid curves. Fig. [2](a) shows the results of 
the xy model with L x = 300/im, L y = 100/im. This 
model describes the in-phase EM emission from IJ Js with 
infinite junction number. The angular distribution of 
the J m has the cylindrical symmetry, and The angular 
dependence of J can be derived from Eq. ([2|) 
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where the 7 = A^/A = 1.36. A^ is the standing wave 
length of E at the boundary, and A is the wave length 
of the far field EM wave. Fig. EJb) shows the radia- 
tion in xz-plane of the sample with L x = 100/im and 
L z = 0.2A C . The thickness of the source is much smaller 
than the wave length. Neglecting the r' in Eq. ([2]), the 
Poynting vector is proportional to sin 2 the dipole angu- 
lar distribution. It is the small size of effective radiation 
section that generates the dipole angular distribution. 



By using a 2D model, Koyama et al [9]. simulated 
the angular distribution of the radiation energy. Their 
numerical results are also list in Figs. [2](a) and (b) by 
open circles. It is very interesting that our analytical 
results agrees quite well with the simulation ones. 

It was recently observed that small samples can be 
used to reduce the heating effect, therefore, the small 
size approach has been generally adopted. When the EM 
radiation from the IJJs is approximately monochromatic 
in THz band, the total power is 
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When the change of the cavity resonance is not obvious, 
(j) x (uj) and (j>t(oj) is constant. Thus Sl^uo 2 / P almost re- 
mains unchanged, where S e ff is the area of the effective 
section. The total power is proportional to the square of 
the number of layers, consistent with the experimental 
observations [4]. 

The power of radiation is also simulated with the uses 
of the sine-Gordon equations and Eq. ([5]). Since the 
magnetic field of the THz EM wave is much smaller 
than the external one, the boundary condition can be 
set as (j) x > = H',H' = 2\ c DH\y^. For such a bound- 
ary condition, the Poynting vector can not be calculated 
directly [l3[, but Eq. ([5]) provides another way to the 
evaluation of the power. The parameters used here are 
Xab = 0.4/im, A c = 200/im, s = 0.3nm, D = 1.2nm, 
/3 = 0.02, L x = 20/im and L y = 20/nm. The number 
of junctions is Nj = 30. The power densities obtained 
in the previous simulations [H, 0] are independent of the 
effective section area. Lin et alfl] estimated the power 
density to be around 400 W/ cm 2 , and Tachiki et al.\^ 
about 3000 W/cm 2 . In sharp contrast with the previous 
reports, Eq. ([5]) shows that the power density is pro- 
portional to S e ff. Our effective section area is set to 
0.9/im 2 , close to the experimental value@. Fig. [3] shows 
the external current dependence of the power intensity 
for the external magnetic field H ext = 2T along the y 
axis. The blue curve denotes the power averaged from 
13.67 ns to 27.34 ns, and the red one from 27.34 ns to 
41.01 ns. Because of the appropriate boundary condi- 
tions we employed, our result for the power density is 
much smaller than previous numerical ones reported in 
the literature 0, Q. The maximum output for the power 
density is about 6 W/cm 2 , which is just in the same or- 
der of magnitude as the experimental result observed by 
Bae et al. p. 

Both Figs. [3ja) and (b) have been divided into two dif- 
ferent regimes, marked with different background colors. 
When J ext < 0.47J C (in the purple regime), the vortices 
are randomly distributed. In general, the power roughly 
increases with J ex t. And if we look into details, it could 
be found that the blue line and red line do not coincide, 
due to an unstable power caused by random flux flow. 
When J ext > 0.47J C , the vortex structure is transformed 
into the in-phase rectangular lattice, a stable configura- 
tion. These configurations are discussed analytically by 
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FIG. 3: (Color online). Radiation power as a function of 
the external current. The blue curve represents the radiation 
power averaged from 13.67 ns to 27.34 ns, the red curve the 
radiation power averaged from 27.34 ns to 41.01 ns. (a) In the 
pink regime, the vortices are distributed almost randomly. In 
the purple regime, the moving Josephson vortices form the 
in-phase rectangular lattice, (b) is an enlarged view of (a). 



Koshelev[14j. In this situation, the distinct power radi- 
ation appears. In addition, the output power fluctuates 
within the pink regime. This phenomenon is closely re- 
lated to the radiation mechanism, and will be discussed 
below. 

It is interesting to note that only a small part of the in- 
put power, not more than 5%, is able to emit out as EM- 
wave in all previous experiments [4|, [6[ . Eq. ([5]) can give 
a reasonable explanation to theses experimental findings. 
When L y Nj is small, the interior surfaces reflect a wave 
of low frequencies. The inner EM- wave, generated by AC 
Josephson effect, bounces back and forth within BSCCO. 
Like a cavity resonator, where the EM field forms stand- 
ing waves, the EM waves with resonance frequencies can 
also form in BSCCO. Every discrete peak in Fig. [3] pro- 
vides an evidence of the standing waves. Fig. [31(b) is an 
enlarged view of Fig. [3](a). It is clear that the distances 
between neighboring peaks are almost the same. These 
peaks just correspond to the characteristic frequencies 
with different resonance modes. 

The different flux-line lattice structures results in the 
different V-I curves and frequency characteristics [ll|. 
To understand the resonance mechanism, we analyze the 
frequency spectra with different external currents. As 
plotted in the Fig. |H[a) and (b), the frequency spectra 
of random lattice is much different from that of the rect- 
angular lattice. 

Fig. |H[a) shows the frequency spectra with external 
currents J ext = 0.21J C , 0.22J C and 0.23J C . Under these 
conditions, the output power is less than 1% of the maxi- 
mum. With low J ex t, the spectra are not monochromatic 
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FIG. 4: (Color online), (a) The frequency spectra of the 
power. The external current are 0.21J c (pink), 0. 2 2 Jc( yellow) 
and 0. 23 J c (blue). Arrows denote the major frequencies. The 
inset shows the enlarged view, (b) Frequency spectra of the 
power with high J ex t- 



and lots of peaks appear. The envelope lines of these 
peaks have characteristic frequencies (black arrow) and 
higher harmonics (red arrows). These frequencies are 
generated from the sin(0). Integration of sine-Gordon 
equations gives ((3 fa) + (sin ((/>)) = J extl where ((3 fa) is 
the normal current and (sin ((/>)) is the super current. As 
i>t) is not zero, the base frequency of sin(0) is derived 
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The frequency, based on the external current, is of great 
importance to tune the supercurrent frequency. Further- 
more, it can be used to control the radiation frequency 
and power. 

Fig. Ufa) shows a discrete cavity resonance and an 
equidistant frequency difference. The phenomenon can 
be clearly seen in the inset, which is an enlarged view. It 
is also interesting to note that in the three lines for differ- 
ent external currents J ext = 0.21 J c (pink), 0.22 Jc(yellow) 
and 0.23J c (blue), the peaks are in the same frequencies. 
It follows that these frequencies are intrinsic and inde- 
pendent of the external parameters, such as the external 
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FIG. 5: (Color online). The intervals of the intrinsic frequen- 
cies. Only one of the parameters(iV s , if, L x , ( and (3) is 
changed in same group of data points. 

current, the magnetic field, and so on. The differences 
(Acj) between two neighboring frequencies could be re- 
garded as the characteristic frequencies. 

In order to clarify these intrinsic frequencies, Aco is cal- 
culated under different J extl H extl length, layer number, 
and coupling C s . The results are divided into 5 groups and 
collected in Fig. [5l The parameters for the same group 
of data are fixed, except that one parameter denoted in 
the legend of Fig. [5] can be tuned. Because the frequen- 
cies are independent of (3 and H, the open circles and 
crosses with different (3 and H coincide. From this fig- 
ure, we can observe that all data point can collapse onto 
a straight line described by Auj = . It follows that 

the intrinsic frequencies is na;i Y^l where N s = Nj + 1 is 
the number of CuO layers and n x is the natural number. 
Through observing the simulation data of the inner EM 
field, it could be found that n x is the number of the stand- 
ing wave nodes. On possible explanation for these phe- 
nomena is the EM resonance between z and x directions. 
The operator (1 — (A^) in Eq. (pQ) could be regarded as 
(1 — d^,) with the dimensionless quantity z' = z e ff/\ c , 
and the effective size z e ff = z\ c / (s + D) The reso- 
nance occurs when uo/k x = u; z /k e ff. 

For high external current, although the higher harmon- 
ics sill exist, the frequency spectra of rectangular lattice 
is almost monochromatic, as shown in Fig. |H[b), close to 



the experiment 1 15| . When the external frequency esti- 
mated by Eq. (|6|) equals to one of the intrinsic frequen- 
cies, the external current would excite a monochromatic 
resonance mode and the radiation would generate the 
maximum power, as denoted by the curves a and b. On 
the other hand, the external current would excite several 
resonance modes, when it is not equal to any intrinsic 
frequency, as shown by the curve c. The pink regime 
of Fig. (j3]) shows these resonance peaks with different 
J ex t- The distance (AJ ext ) between two neighbor peaks 
is 0.0098 J c , consistent with the numerical result 0.0093 J c 
based on intrinsic frequencies. 

For fields parallel to the ab— plane in BSCCO, vortices 
are periodically distributed. The distributions of super- 
current and the phase of the order parameter are also 
periodic. The normal current is no longer uniform, which 
results in periodical perturbation in the inner electric 
field. Therefore the vortices excite the alternative EM 
field. When the vortices are in motion, the alternative 
EM field is also propagating and decaying in BSCCO. In 
order to form a stable alternative EM field, the vortices 
and the EM field must move synchronously. 

By increasing or decreasing the external current, the 
velocity of the flux flow is also changed. When J ext < 
0.47J C , the flux flow would move more slowly than the 
EM field. To achieve synchronies along the x-axis, the 
normal direction of the EM wavefront is not parallel to 
the flow direction. Because the boundary surface and 
the wavefront is not tangential, which weakens the power 
emission. On the other hand, when J ext > 0.47J C , the 
external frequency would possibly resonate with the in- 
trinsic frequencies, which pumps considerable energy into 
the inner alternative EM field. Since the boundary sur- 
face is tangential to the wavefront, the power emits much 
more efficiently. 
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